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Abstract. One of the most effective techniques of experimental mathematics 
is to compute mathematical entities such as integrals, series or limits to high 
precision, then attempt to recognize the resulting numerical values. Recently 
these techniques have been applied with great success to problems in math- 
l | ematical physics. Notable among these applications are the identification of 

("^ some key multi-dimensional integrals that arise in Ising theory, quantum field 

O | theory and in magnetic spin theory. 

•s 



1. Introduction 



One of the most effective techniques of experimental mathematics is to compute 
mathematical entities to high precision, then attempt to recognize the resulting 
numerical values. Techniques for efficiently performing basic arithmetic operations 
and transcendental functions to high precision have been known for several decades, 

r-4- and within the past few years these have been extended to definite integrals, sums 

of infinite series and limits of sequences. Recognition of the resulting numerical 

ly-^ values is typically done by calculating a list of n possible terms on the right-hand 

side of an identity, also to high precision, then applying the PSLQ algorithm [211 111] 
to see if there is a linear relation in this set of n + 1 values. If PSLQ does find a 
credible relation, then by solving this relation for the value in question, one obtains 
£> a formula. These techniques have been described in detail in |14j . |15j . and [9]. 

JJ~j In almost applications of this methodology, both in sophistication and in com- 

r% putation time, the most demanding step is the computation of the key value to 

sufficient precision to permit PSLQ detection. As we will show below, computation 
of some high-dimensional integrals, for instance, often requires several hours on 
a highly parallel computer system. In contrast, applying PSLQ to find a relation 
among, say, 20 candidate terms, each computed to 500-digit precision, usually can 
be done on a single-CPU system in less than a minute. 

In our studies of definite integrals, we have used either Gaussian quadrature 
(in cases where the function is well behaved on a closed interval) or the "tanh-sinh" 
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quadrature scheme due to Takahasi and Mori 29 (in cases where the function has 
an infinite derivative or blow-up singularity at one or both endpoints) . For many in- 
tegrand functions, these schemes exhibit "quadratic" or "exponential" convergence 
- dividing the integration interval in half (or, equivalently, doubling the number of 
evaluation points) approximately doubles the number of correct digits in the result. 
The tanh-sinh scheme is based on the observation, rooted in the Euler-Maclaurin 
summation formula, that for certain bell-shaped integrands (namely those where 
the function and all higher derivatives rapidly approach zero at the endpoints of 
the interval), a simple block- function or trapezoidal approximation to the integral 
is remarkably accurate pg. 180]. This principle is exploited in the tanh-sinh 
scheme by transforming the integral of a given function f(x) on a finite interval 
such as [— 1, 1] to an integral on (— oo, oo), by using the change of variable x = git), 
where g(t) = tanh(7r/2-sinht). The function g(t) has the property that g(x) — > 1 as 
x — > oo and g(x) — > —1 as x — Y — oo, and also that g'(x) and all higher derivatives 
rapidly approach zero for large positive and negative arguments. Thus one can 
write, for h > 0, 

/l ,.00 n 

f{x) Ax = / f(git))g'(t)dt « h £ t^/fo), 
-1 J-oo j=-N 

where the abscissas Xj — g(hj), the weights Wj — g'ihj), and N is chosen large 
enough that terms beyond N (positive or negative) are smaller than the "epsilon" 
of the numeric precision being used. In many cases, even where fix) has an infinite 
derivative or an integrable singularity at one or both endpoints, the transformed in- 
tegrand figit))g'(t) is a smooth bell-shaped function for which the Euler-Maclaurin 
argument applies. In these cases, the error in this approximation (nl) decreases more 
rapidly than any fixed power of h. Full details are given in |12j . 

Both Gaussian quadrature and the tanh-sinh scheme are appropriate for an- 
alytic functions on a finite interval. Functions on a semi-infinite intervals can be 
handled by a simple transformation such as: 



fit) At = [ fit)At+ [ 
o Jo Jo 



/(lA)dt 
t 2 



Oscillatory integrands such as J Q (l/xsina;) p da: can be efficiently computed by 
applying a clever technique recently introduced by Ooura and Mori |26j . Let x — 
git) = Alt /il — cxp(— 27rsinhi)). Then in the case of p — 2, for instance, 

sm g(t)\ , 

■ g it) At 



rm * - / 



git) 

-^ fsmgihk)\ 2 , 



Now note that if one chooses M — ir/h, then for large k, the g{hk) values are all 
very close to fc7r, so the sin(g(/iA:)) values are all very close to zero. Thus the sum 
can be truncated after a modest number of terms, as in tanh-sinh quadrature. In 
practice, this scheme is very effective for oscillatory integrands such as this. 
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In the next four sections we consider Ising integrals, Bessel moment integrals, 
'box' integrals, and hyperbolic volumes arising from quantum field theory respec- 
tively. We then conclude with a description of very recent work on multidimensional 
sums: Euler sums and MZVs. 

2. Ising integrals 

In a recent study, Bailey, Borwein and Richard Crandall applied tanh-sinh 
quadrature, implemented using the ARPREC package, to study the following classes 
of integrals [8]. The D n integrals arise in the Ising theory of mathematical physics, 
and the C n have tight connections to quantum field theory. 

4 f°° f x 1 du\ du„ 



D„ 



n - JO JO fx^ n („, i 1 /„, ~iA u l u r, 



(E^=i(«i + V«i))' 

iu<j y-Ui+Uj J dui du r , 

(E-=i(^ + V 



En = 2 [■■■[( n hs ^^-*- 



3) 

1 



2 U\ u r , 



J ° \l< 3 <k<n Uk 



where (in the last line) u^ — IIi=i ^i- 

Needless to say, evaluating these n-dimensional integrals to high precision 
presents a daunting computational challenge. Fortunately, in the first case, we were 
able to show that the C n integrals can be written as one-dimensional integrals: 



C n = -r pK(p)d P , 

n\ Jo 



where Kq is the modified Bessel function jl]. After computing C n to 1000-digit 
accuracy for various n, we were able to identify the first few instances of C n in 
terms of well-known constants, e.g., 

Cs = L - 3(2) = g ((^TTF - WTw) 
Ci = i^ c(3) ' 

where £ denotes the Riemann zeta function. When we computed C n for fairly large 
n, for instance 

Ci 024 = 0.63047350337438679612204019271087890435458707871273234 

we found that these values rather quickly approached a limit. By using the new edi- 



tion of the Inverse Symbolic Calculator, available at http://ddrive.cs.dal.ca/~isc 
this numerical value can be identified as 

lim C n = 2e- 27 , 

n— >oo 

where 7 is Euler's constant. We later were able to prove this fact — this is merely 
the first term of an asymptotic expansion — and thus showed that the C n integrals 
are fundamental in this context [5]. 
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The integrals D n and E n are much more difficult to evaluate, since they are 
not reducible to one-dimensional integrals (as far as we can tell), but with certain 
symmetry transformations and symbolic integration we were able to reduce the 
dimension in each case by one or two. In the case of D§ and E§ , the resulting 3-D 
integrals are extremely complicated, but we were nonetheless able to numerically 
evaluate these to at least 240-digit precision on a highly parallel computer system. 
In this way, we produced the following evaluations, all of which except the last we 
subsequently were able to prove: 

D 2 = 1/3 

D 3 = 8 + 4tt 2 /3-27L_ 3 (2) 

At = 4tt 2 /9 - 1/6 - 7C(3)/2 

E 2 = 6 - 8 log 2 

E 3 = 10 - 2tt 2 - 8 log 2 + 32 log 2 2 

Ei = 22-82<C(3) -241og2 + 1761og 2 2-256(log 3 2)/3 + 167r 2 log2-227r 2 /3 

E 5 = 42-1984Li 4 (l/2) + 1897r 4 /10-74C(3)-1272C(3)log2 + 407r 2 log 2 2 
-62tt 2 /3 + 40(tt 2 log 2)/3 + 88 log 4 2 + 464 log 2 2-40 log 2, 

where Li denotes the polylogarithm function. In the case of D 2 , D3 and D4, these 
are confirmations of known results. We tried but failed to recognize D 5 in terms of 
similar constants (the 500-digit numerical value is available if anyone wishes to try). 
The conjectured identity shown here for E 5 was confirmed to 240-digit accuracy, 
which is 180 digits beyond the level that could reasonably be ascribed to numerical 
round-off error; thus we are quite confident in this result even though we do not 
have a formal proof. 

In a follow-on study [6], we examined the following generalization of the C n 
integrals: 

• , , 1 dwi Au n 



" !7 " Jo (E"=i(^ + i/^)) fc+1 Ul ""' 

Here we made the initially surprising discovery — now proven in |17j and in outline 
much earlier |13) — that there are linear relations in each of the rows of this array 
(considered as a doubly-infinite rectangular matrix), e.g., 

= C 3 . - 84C 3 , 2 + 216C 3 ,4 

= 2C 3 ,i - 69C 3 , 3 + 135C 3 , 5 

= C za - 24C 3 , 4 + 40C 3 , 6 

= 32C* 3 , 3 - 630C 3 , 5 + 945C 3 , 7 

= 125C 3 . 4 - 2172C 3 , 6 + 3024C 3 , 8 . 

3. Bessel moment integrals 

In a more recent study of Bessel moment integrals, co-authored with Larry 
Glasser [7], the first three authors were able to analytically recognize many of 
the C n ,k constants in the earlier study — because, remarkably, these same integrals 
appear naturally in quantum field theory (for odd k). We also discovered, and then 
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proved with considerable effort, that with c n ^ normalized by C n> k — 2™ c n ^/(n\ kl), 
we have 



C3,0 



C3,2 



C4,0 



C4,2 



3r 6 (l/3) VSit 3 



32tt2 2 / 3 


288 \ 


4 oo /2n\ 4 



4 Z-^i 4477 

n=0 



8 "^ { 

1/2,1/2,1/2 
2,2 



-4-^3 



1/2,1/2,1/2 
1,1 



1/2,1/2,1/2,1/2 



7T 

64 



4 4 F; 



4^3 



1/2,1/2,1/2,1/2 
1,1,1 



-3 4 F 3 



1/2,1/2,1/2,1/2 
2,1,1 



1 



3tt^ 
16 ' 



where p F q denotes the generalized hypergeometric function [1 . The corresponding 
values for small odd second indices are 03,1 = 3L_ 3 (2)/4, 03,3 = L_ 3 (2) — 2/3, 04,1 = 
7C(3)/8 and c 4 , 3 = 7C(3)/32 - 3/16. 

Integrals in the Bessel moment study were quite challenging to evaluate nu- 
merically. As one example, we sought to numerically verify the following identity 
that we had derived analytically: 

it r' 2 r' 2 K(sinfl)K(sin0) 

C5,o = 7; I / / dOdcj), 

-tt/2 J-1T/2 , / cos 2 Q cos 2 A , 4 sin 2/0 , ja 



where K denotes the elliptic integral of the first kind pQ . Note that this function 
has blow-up singularities on all four sides of the region of integration, with par- 
ticularly troublesome singularities at (n/2,— it/2) and (— 7r/2, tt/2) (see Figure 1). 
Nonetheless, after making some minor substitutions, we were able to evaluate (and 
confirm) this integral to 120-digit accuracy (using 240-digit working precision) in a 
run of 43 minutes on 1024 cores of the "Franklin" system at LBNL. 

In a separate study, the first two authors studied correlation integrals for the 
Heisenberg spin-1/2 antiferromagnet, as given by Boos and Korepin, for a length-n 
spin chain [24, eqn. 2.2]: 



P{n) 



OO poo 



7r n(n+l)/2 

(2ni) n 
x T(xi - i/2, x 2 - i/2, 



— OO J — OO 



U(xi-i/2,x 2 -i/2,-- 

-CO 

, x n - i/2) dxi dx 2 • • • dx n 



i/2) 



where 
U(x!-i/2,x 2 -i/2, 

T(xi - i/2,X2 - i/2, 



j X% 



i/2) 
i/2) 



Yii 



<k<j<n 



sinh[7r(iEj — Xk)] 



Ui< 3 <n inc osh n (n Xj ) 
ni<,-< n (^--i/2) J '- 1 (^- + i/2) ft - J ' 



III 



<fc<i< 



n ( X 3 



Xk 



They computed numerical values for these n-fold integrals to as great a pre- 
cision as we could, then attempted to recognize them using pslq. They found 
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Figure 1. Plot of c$ t o integrand function 



the following, which confirm some earlier results obtained by others using physical 
symmetry methods: 



(1) = 


1 
2 


(2) = 


= log2 

3 3 B 


(3) = 


= 4-log2 + 


(4) = 


= i-21og2 + 
5 


(5) = 


1 io, „ 

= log 2 

6 3 B 




1225 w N , 



|«3) 

f«3)->).o g 2-^(3)-g«5) + 



8,-) 



^ C (3)-^C(3)log2-^C 2 (3)-^ 

24 SV I 2 SV ! 5 lg S V ! 1Q2 



24 
C(5) 



C(5)log2 



425 



12125 



6 C(5)log2-^ C (3)C(5)- 256 

11515 C(7)log2 + ^ C (3)C(7) 



C 2 (5) 



6223 
~256 



C(7) 



64 



512 



These computations underscore the rapidly increasing cost of computing inte- 
grals in higher dimensions. Precision levels, processor counts and run times are 
shown in Table 1. 
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n 


Digits 


Processors 


Run Time 


2 


120 


1 


10 sec. 


3 


120 


8 


55 min. 


4 


60 


64 


27 min. 


5 


30 


256 


39 min. 


6 


6 


256 


59 hrs. 



Table 1. Run times and precision levels for spin integral calculations 



4. Box integrals 

Let us define box integrals for dimension n as 



B n (s) 
A n (s) 



s/2 



dri ■ ■ ■ dr„ 



s/2 



((ri - qi) 2 -\ h (r„ - q n ) 2 ) dri • • • dr„ dq ± ■ ■ ■ dq n . 



(2) 
(3) 



As explained in previous treatments [4], [5], these integrals have several physical 
interpretations: 

(1) B n (l) is the expected distance of a random point from the origin (or from 
any fixed vertex) of the n-cube. 

A n (l) is the expected distance between two random points of the n-cube. 
B n (—n + 2) is the expected electrostatic potential in an n-cube whose 
origin has a unit charge. Such statements presume that electrostatic po- 
tential in n dimensions is V(r) = 1/r" -2 , and say logr for n = 2; in other 
words, the negative powers of r can also have physical meaning. 

(4) A„(— n + 2) is the expected electrostatic energy between two points in a 
uniform cube of charged "jellium." 

(5) Recently integrals of this type have arisen in neuroscience e.g., the average 
distance between synapses in a mouse brain. 

Note that the definitions show immediately that both A n (2m) and B n (2m) are 
rational when m, n are natural numbers. A pivotal, original treatment on box 
integrals is the 1976 work of Anderssen et al [2]. There have been interesting 
modern treatments of the B n and related integrals, as in |10j . |14l pg. 208], |32j . 
and |30j . Related material may also be found in |23L 131 j . 

Like the Ising integrals, some of these n-dimensional integrals are reducible to 
1-dimension integrals. For instance, we found that 



A 8 (-l) = 



(-l + e~" + % /7ruerf(w)) 3 



du. 



After calculating a 400-digit numerical value for this constant, we were able to 
recognize it as 



A 3 (-l) = — (6 + 6\/2-12\/3-107r + 301og(l 
15 V 



/ 2) + 301og(2 + \/3) 



A selection of results that we have found are shown in Tables 2, 3, 4 and 5. 
Here G denotes Catalan's constant, namely G :— ^ n>0 {—^-) n / {^n + l) 2 , 9 — 
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n 


s 


B n (s) 


any 


even s > 


rational, e.g.: i?2(2) = 2/3 


1 


,s^-l 


l 

s+l 


2 


-4 


4 8 


2 


-3 


-v/2 


2 


-1 


21og(l + y2) 


2 


1 


y2+\\og{l + V2) 


2 


3 


^x/2+|jlog(H-%/2) 


2 


s^-2 


2+s 2i ! V2' 2' 2' ■"■/ 


3 


-5 


_ 6 ^3 - i^TI" 


3 


-4 


— | \[2 arctan -4= 


3 


-2 


-3G + |tt Iog(l + \/2) + 3Ti 2 (3 - 2\/2) 


3 


-1 


-i7r+|log(2 + \/3) 


3 


1 


£73-^+1 log (2 + ^3) 


3 


3 


1^-^+^^(2 + 73) 



Table 2. Recent evaluations of Box integrals 



n 


s 


B n {s) 


4 
4 
4 
4 

4 


-5 
-3 

-2 

-1 

1 


— \/8 arctan 1 -4= J 

4G-12Ti 2 (3-272) 

7T log (2 + V3) - 2 G - £ 

2 log 3 - § G + 2 Ti 2 (3 - 2 72) - 78 arctan (^) 

I - f§ + io Ti 2 (3-2^2) +log3- ^ arctan(^) 


5 

5 
5 

5 


-3 

-2 
-1 

1 


^ G - 10 log (2 - V3) $ - | tt 2 

-10 Cl 2 (| + | 7r) + 10 Cl 2 (| e - | 7r) + f Cl 2 (^ + | 7r) 

+ fCl2(^+! 7 r)-fCl2(^ + | 7 r)-fCl2(e + f.) 
| J B 5 (-6)-|B 5 (-4) + |7rlog3 + 10Ti 2 (i) -10G 

-iMG+f^log(2- v / 3) + i^ 

+5 log ( 1+ /^) 2 ^ arctan ( ^) 

_f Cl 2 (0+§7r) + f Cl 2 (0+§7r) + f Cl 2 (<?+lH 

- g G + I 9 log (2 - 73) + 3^o ^ 2 + | V5 
+ f ^g (^)~|V3 arctan (-^) 

+ 7 Cl2 (l < , + l 7r )_7 Cl2( l < ,_l 7r) _> Cl2(e+ l 7r) 

_i4ci 2 (^+|7r) + ^Cl 2 (^+|7r) + i|Cl 2 (0+^7r) 



Table 
arctan 



3. Recent evaluations of Box integrals, continued; here i 

16-3\/l 5 N 

n 
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n 


s 


A„(s) 


2 


-5 


^ + ^/2 

3 + 9 VZ 


2 


-1 


f-|v^ + 41og(l + >/2) 


2 


1 


^ + i^+ilog(l + y2) 


3 


-7 


4 16V2 | 2V3 , 7T 

5 15 + 5 + 15 


3 


-2 


2tt - 12 G + 12 Ti 2 (3 - 2^2) + 6tt log (l + \/2) 
+2 log 2- § log 3 -8^2 arctan f-L) 


3 


-1 


| fTr+fv^ |v / 3 + 21og(l + ^) + 121og( 1 ^) 41og(2 + V3) 


3 


1 


"f 1^+21^ |v / 3 + 2 1og(l + y2)+8 1og( 1 ^f) 


3 


3 


-i^-:E55*+ii>v^+^+|stog(i + v^ + itog(±±^) 



Table 4. Recent evaluations of Box integrals, continued 



A„(s) 



-1 

1 



128 

15 
8 
5 



-IV2 



f 7T - 8 log (1 + v/2) - 32 log (1 + v/3) + 16 log2 
3-32\Z2arctanf-^") - 96 Ti 2 (3 - 2 >/2 
f 7rlog2+^ tt log (1 + ^3 



3 



152 8_ 

315 15 
1 32 

"•" 5 

23 
135 



+ f \/3- 32 \/2 arctan (^ 
i|7ry3-|^lo g 2+f 7rlog(l 
+ | V2 arctan (2 y/2) + § log 3 - 4 tt log (\/2 - l) 
+8Ti 2 (3 - 2 s/2) - f G - § log 2 
log2+|log3+^|y2-i|y3 
G + 8Ti 2 (3 - 2 ^2) - § 



20 log 3 
f 32 G 



I- 



TT 

log(l 



ill 
a/3) 



log (1 + 7 
v^ arctan (^2/4) 



16 - 52 1082+^1083+^^+^^3 



315 



7T 



105 



log (1 + x/2) 



Y^ \[2 arctan 



420 
1 104 
"^ 105 
4 /o 
~~ 15 U 



630 

log (1 + \/3) 



105 



Ti 2 (3 - 2 V2 



1279 G __4_,^_ 7r 2_J49_, 3239 /o 
567 ""* 189 ~ 315 3465 ~ 62370 v 

ilog(l + ^) + ilog(2 + y3) 

^log(2-^)^ 



568 /o _ 380 
3465 v ° 6237 

' log (±Vy 



_380_ ./c 
6237 v u 



64 

18 



-§§§ log 3 -r 54 iug v ± -r v */ t gg iug v * -r v «V -r m i"i 

"I V2 arctan (^) - £ ^3 arctan (^) + ±f log (2 
v / 2) + ^Cl 2 (^+|7r)-^Cl 2 (^ 

04 m f'/i , 1 _\ 208 m f'/i , 4 _\ 



-|Ti 2 (3- 



63 , v „ , 

- 2\/2) -r -gg- >^i 2 ^g v -r 3 /i ; - , : , 

+ I jci 2 (fl + H + Jic^+ 6 H 

Table 5. Recent evaluations of Box integrals, continued 



^ 



arctan((16 — 3vl5)/ii), CI denotes Clausen's function, 

am = E=P. 
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and Ti denotes Lewin's inverse-tan function, 

Ti 2 (x) = £(-1)' 



x 2n+l 



(2n+l) 2 

n>0 V ; 



5. Clausen functions and hyperbolic volumes 

In an unpublished 1998 study |16j two of the present authors (Borwein and 
Broadhurst) identified 998 closed hyperbolic 3-manifolds whose volumes are ratio- 
nally related to Dedekind zeta values, with coprime integers a and b giving 

(2) ?vo>(M) - t-D?»W2) 



b y ' (2tt) 2 »- 4 2C(2) 

for a manifold Ai whose invariant trace field K has a single complex place, discrim- 
inant D, degree n, and Dedekind zeta value C,k (2). While the existence of integers 
a, b can be established, via algebraic if-theory as in |35j . for the most part it was 
and is not possible to specify the rational a/b other than empirically |35j . 
The simplest identity implicit in (pi devolves to 

(3) 3Cl 2 (a)- 3Cl 2 (2a) + Cl 2 (3a) = ^L_ 7 (2), 

with a = 2arctan(v / 7), as is recorded in [141 p. 91]. Here L^{2) := Yl n >o (?) / n2 
is the primitive Dirichlet L-series modulo 7 evaluated at 2 where (?) is the Legendre 
symbol. This was rewritten in equivalent and more self-contained form as 



24 /' 7r/ ' 2 

(4) 3? L bs 



tan t + v7 
tan t — v7 



dt = L_ 7 (2) 



in [9l p. 61] — and elsewhere. 

Note that the integrand function of |4j) has a nasty singularity at arctan(\/7) 
(see Figure p|. However, we were able to numerically evaluate this integral to 
20,000-digit accuracy, by splitting the integral into two parts, namely on the inter- 
vals [n/3, arctan(v / 7)] and [arctan(v / 7),7r/2]. Note that tanh-sinh quadrature can 
be used on each part, since it can readily handle blow-up singularities at one or both 
endpoints of the interval of integration. This run required 46 minutes on 1024 CPUs 
of the Virginia Tech Apple cluster. The right-hand side was also evaluated, using 
Mathematica, to 20,000-digit precision. The two values agreed to 19, 995 digits [9, 
pg. 61]. Alternative representations of the integral in Q are given in |20j . 

We shall now provide a proof of Eqn. (pi) and hence of Eqn. Q . Actually, an 
equivalent (if not obviously so) form of identity (I3J) , namely 



(5) Cq(^t)(2) - ^{A(cot^)+^(coty)+A(coty)} 

- ; {2A(V7) + A(V7 + 2\/3) + A(V7 - 2^)} 



7V7 



with the notation 



f x 1 4 

A(x) ■- / ^-^ log^-p-^dt = Cl 2 (2arccota;) 
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FIGURE 2. Plot of integrand function in Q 



is already established in |33j . The first equality in (pi) can be written as 



(6) 



On noting that 



CQ(y=7)(2) 



L_ 7 (2) = C(2)L_ r (2). 



cot arg 



y/l + % 

C0targ ^7T 

(1+ »>/?)(! T*V3) 



4\/2 






C\ 2 {6) = lmJ2 



n>l 



we can translate the remaining, highly non-trivial, part of ([5]) to 
2A(Vf) + A(V7 + 2-\/3) + A(V7 - 2\/3) 



Z— • TJ^ 



»V V 2V2 



1 ( (i + iV7)(i - % ^iy 2 " 



4V2 



m 4-; « 2 v 4V2 
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Now we use 

2 



(l + n/7)(l + V3) V 2W3 

(l + tV7)(l-i>/3) V _ _ 2W3 



4V2 ' ^ 



/\/7 + «V" T /3 + i\/7 



Im l 2VrJ = M - Y~J = _Im (-/*) n fOT ™ = 0,1,2,..., 

where [i := (—3 + iyf)/4 has absolute value 1 and arg/i = a = 2arctan(y / 7), to 
write the latter equality as 

(7) 2A{V7) + A(V7 + 2V3) + A{V7 - 2V3) 

^ fi"(e 2 ™/ 3 + e - 2 ««/3 _ 2(-l) n ) 



n 

n>l 



2 



V^-^ n 2 ^— ' n 2 3 ^— ' n 2 

Nl>l n>l n>l 

- Cl 2 (a) - Cl 2 (2a) + ^Cl 2 (3a), 

where we have applied the following two standard identities 
l^i 2 " _ ^ x n y, {-l) n x n 

2 2-^1 n 2 2-^1 n 2 Z^ n 2 
n>l n>l n>l 

I ™3n _n JlirinlZ~n — 27rin/3~,n 

1 v -v X v -v X y -v C ' X y -s c ' X 

3 ^-' Ti 2 2^/ n 2 2-^i n 2 Z^ n 2 
n>l n>l n>l n>l 

for the dilogarithm function. It remains to substitute our finding (Fn) into (|5| and 
(|6| to finish a proof of identity ([3]) . 

The equivalent identity Q can be obtained by some reasonably straightforward 
but tedious manipulation of the Clausen integral representation 



(8) 



Cl 2 (6>) =- / log 1 2 sin cr | do- 
Jo 



for < 9 < 2ir, and an appropriate change of variables. 

As Don Zagier points out in |33j 

"we observe that the values of A(x) at algebraic arguments satisfy 
many non-trivial linear relations over the rational numbers; I 
know of no direct proof, for instance, of the equality of the right- 
hand sides of Eqns. (5) and (6)." 

Zagier's Eqns. (5) and (6) are our identity ^. Another result in |33j . Theorem 3, 
implies the identity 



<9) ^> <2) " ^M^+Kt?)^) 
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which may be thought of as complimentary to Eqn. ^ (see pg. 300 in [33] for 
details). Since 

a(-/=) = Cl 2 (2 arctan Vf) 
-A (-7=) = Cl 2 ( 2 arctan^ 

^(-7=) = Cl 2 (2 arctan — 

and 

/7 fl 

2 arctan v7 = a, 2 arctan — = —2a + 27r, 2 arctan — — 3a — 2tt, 

3 5 

identity p| follows from |9| immediately. Thus paper |33j contains two different 
proofs of ([3])! 

Let us clarify the current status and somewhat-complicated history of various 
of the discoveries in |16j . Until recently the authors of |16j after discussion with 
Zagier believed ([5]) to be unproven. It was only when Zudilin spent time with Don 
Zagier in 2008 that he remembered his equivalent pre-dilogarithm (see |34L 135] ) 
result in [33] . Two of the present authors (Borwein and Broadhurst) |16j wrote 

" While the existence of such relations is understood 33l 1341 135] , 
their precise forms appear to be unpredictable, thus far, by de- 
ductive mathematics. They are therefore ripe for the application 
of experimental mathematics." 

The great bulk of the results recorded in [16 remain unproven. They were dis- 
covered by intensive physically and mathematically motivated computation, using 
SnapPea, Pari-GP, Maple, and other tools. 

Indeed, the cases 

D = -8,-11,-15,-20,-24 

are challenging enough! These five respectively yield the following conjectured 
identities — each of which is open. First 

(10) 27C1 2 (0 2 ) - 9C1 2 (20 2 ) + Cl 2 (30 2 ) = 8C1 2 (^- J + 8C1 2 '*'" 
with 82 '■= 2 arctan y/2. Secondly, 

(11) 15Cl 2 (0ii)-lOCl2(20 u ) + Cl 2 (50ii) = ll^^A^ci,^ 

where 8n :— 2 arctan yll and (^-) is the Jacobi (or Legendre) symbol for the 
Dirichlet character. Thirdly, 

(12) 24Cl 2 (0 5j3 ) - 12C1 2 (20 5 . 3 ) - 8C1 2 (30 5 , 3 ) + 6C1 2 (40 5:3 ) 



1 *im^m 
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with #5^ := 2 arctan yJh/Z. Fourthly 

(13) 36C1 2 (# 5 ) - 3OC1 2 (20 5 ) + 4C1 2 (30 5 ) + 3C1 2 (4# 5 ) 

= 20 Yj C1 2 

fcG{l,3,7,9} 

with 9 5 := 2 arctan y5- Finally, 

(14) 6OC1 2 (0 3 , 2 ) - 18C1 2 (2#3, 2 ) - 4C1 2 (30 3 , 2 ) + 3C1 2 (40 3 , 2 ) 

= 24 J2 Cl 2 

fce{i,5,7,n} 



2irk 
~20~ 



2irk 
~24 



with 6*3 >2 := 2 arctan a/3/2. So, for the fifth time, we have a relation that is as 
easy to check numerically as it appears hard to derive. Needless to say, it would 
be interesting to check whether Zagier's 1986 theorems in |33) work for all such 
small values of D; Theorem 2 in 33 looks sufficiently powerful for this task, while 
Theorem 3 therein depends critically on a delicate geometric construction and might 
be of use for D = —11, —15, —20. Moreover, is there a more transparent method to 
deduce identity ^ as well as |lOJHfl4|)? 

As another example of the ubiquity of Clausen values, we complete this section 
with the most difficult integral evaluation required in [5J: 



(15) /C, := / 4 ; TCSeC{x) dx 



3 n, 6 M 3 m ( 6 M 18 , 

3 ) - n C H3 " 6 ) - H C H3 + 6 ) + IT' 

^ci 2 (e- 2 -f\- -ci 2 (* + *)- ^ci 2 (e-l 

11 V 3 / 11 V 6/ 11 V 6 



(»-5) 



20-|)log(2-v / 3) 



Here 



/ 16-3V15 

:= arctan 

11 



It may well be that this closed form ( 15 ) for K,\ can be further simplified 



6. Relations between MZVs and Euler sums 

We conclude with an application of experimental mathematics to discover re- 
lations between multiple zeta values (MZVs) of the form 

1 



C(si,s 2 ,...,s fc ) = Y 



si Su 

n, ... n h 

ni>n 2 >...>n k >0 L K 



with weight w — y^-_ 1 Si and depth k and Euler sums of the more general form 

E 



jn i 



7li 
ni>n 2 >...>n k >0 1 
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with signs a = ±1. Both types of sum occur in evaluations of Fcynman diagrams 
in quantum field theory [181 119) as mentioned in [14] . These sums are described 
in some mathematical detail in [15, Chapter 3]. 

First we recall the first Broadhurst-Kreimer conjectures (see |18j and also |15j ) 
for the enumeration of primitive MZVs and Euler sums of a given weight and depth. 
Let E n< k be the number of independent Euler sums at weight n > 2 and depth k 
that cannot be reduced to primitive Euler sums of lesser depth and their products. 
It is conjectured that 



nTT fl - x n v k ) B "- k = 1 — 

We emphasise that, since the irrationality of odd values of depth-one MZVs 
(i.e., Ricmann's £) is not settled, such dimensionality conjectures are necessarily 
experimental. Now let D n ^ be the number of independent MZVs at weight n > 2 
and depth k that cannot be reduced to primitive MZVs of lesser depth and their 
products. Thus we believe that Z?i2,4 = 1, since there is no known relationship 
between the depth-4 sum C(6,4, 1, 1) = J2j>k>i>m ^-/{j^kHm) and MZVs of lesser 
depth or their products. It is conjectured that 



nna-^T-' = i- 



x A y x y (1 — y ) 



l-x 2 (1 - x A )(l - x 6 )' 

»>2fc>0 v /v ; 

The final Broadhurst -Krcimcr conjecture concerns the existence of relations 
between MZVs and Euler sums of lesser depth. The now proven relation |19| 

C(6,4,l,l) = ^ C (9,3) + ^C(9,3) + 3C(2)C(7,3) + ^C(5)C(7) 

1 Q1Q1 7 

+ ^C 4 (3) - W C(3)C(9) + 2C(2)C 2 (5) + 10C(2)C(3)C(7) 



-C 2 (2) 



^C(5,3)4C(3)C(5)-|C(2)C 2 (3)-^|C 4 (2) 



shows that the depth-4 MZV on the left can be expressed in terms of Euler sums of 
lesser depth and their products. In fact, it suffices to include the alternating double 
sum C(9, 3) = J2j>k>o(~ l) J ' +fe /(_7 9 fc 3 ), where a bar above an argument of C serves 
to indicate an alternating sign. In the language of |18l 119] this is a "pushdown" , 
at weight 12, of an MZV of depth 4 to an Euler sum of depth 2. Let M n> k be 
the number of primitive Euler sums of weight n > 2 and depth k whose products 
furnish a basis for all MZVs. It is conjectured that |18j 

nn^-^r^ = i x3y 



i 

n>2fc>0 



Then by comparison of the output -D2i,3 == 6, D21.5 — 9, i?2i,7 — 1 of ( 16 1 with the 



output M21 3 = 9, M21 5 = 7 of (16) we conclude that at weight 21, for example, 
three pushdowns are expected from depth 5 to depth 3 and one from depth 7 to 
depth 5. 

By massive use of the computer algebra language FORM, to implement the 
shuffle algebras of MZVs and Euler sums, the authors of |19j were recently able to 
reduce all Euler sums with weight w < 12 and all MZVs with w < 22 to concrete 
bases whose sizes are in precise agreement with conjectures ( 16|16 ). Moreover, 
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further support to these conjectures came by studying even greater weights, w < 30, 
using modular arithmetic. However, such algebraic methods were insufficient to 
investigate pushdown at weight 21. Instead the authors resorted to a combination 
of the PSLQ methods reported in [IT] with the lll algorithm [25 of Pari-GP 27J, 
finding empirical forms for precisely the expected numbers of pushdowns at all 
weights w < 21. Most notable of these is the pushdown from depth 7 to depth 5, 
at weight 21, in the empirical form 

QO/? 

C(6,2,3,3,5,l,l) = --r-C(3, 6, 3, 6, 3) + {depth- 5 MZV products} 

81 

where the remaining 150 terms are formed by MZVs with depth no greater than 5, 
and their products. 

It is proven, by exhaustion, in |19j that the shuffle algebras do not allow the 



sum £(6,2,3,3,5, 1, 1) in equation (16 1 to be reduced to MZVs of depth less than 
7. It is also proven that all other MZVs of weight 21 and depth 7 are reducible to 
£(6, 2, 3, 3, 5, 1, 1) and MZVs of depth less than 7. Yet it appears to be far beyond 
the limits of current algebraic methods to prove that inclusion of the rather striking 
depth-5 alternating sum 



C(3,6,3,6,3) = x 



^— ^ (jk 2 lm 2 n) 3 ' 

j>k>l>m>n>0 w ' 



with the rather simple coefficient —326/81, leaves the remainder reducible to MZVs 
of depth no greater than five. 

Thus we are left with a notable empirical validation of a pushdown conjecture 
relevant to quantum field theory, crying out for elucidation. 



7. Conclusion 

We have presented here a brief survey of the rapidly expanding applications 
of experimental mathematics (in particular, the application of high-precision arith- 
metic) in mathematical physics. It is worth noting that all but the penultimate of 
these examples have arisen in the past five to ten years. Efforts to analyze inte- 
grals that arise in mathematical physics have underscored the need for significantly 
faster schemes to produce high-precision values of 2-D, 3-D and higher-dimensional 
integrals. Along this line, the "sparse grid" methodology has some promise |28M36j . 

Current research is aimed at evaluating such techniques for high-precision ap- 
plications. To illustrate the difficulty, we leave as a challenge to the reader the 
computation of the triple integral 

v / f(u,v,w)-2dudvdw = 1.1871875..., 
c 

where C := [0, 1/2] 3 and 

f(u,v,w) := cos 2 ((v + w)ir) + cos 2 ((u — v)ir) + cos 2 ((u + w)w) 
+ cos 2 (vir) + cos 2 (u7r) + cos 2 (wit) 

to, say, 32 decimal digit accuracy. 
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